Fermionization Transform for Certain Higher-Dimensional Quantum Spin Models 
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Proposed is a generalization of Jordan- Wigner transform that allows to exactly fermionize a large family 
of quantum spin Hamiltonians in dimensions higher than one. The key new steps are to enlarge the Hilbert 
space of the original model by adding to it a collection of stand-alone free spins and to use a combination 
of these auxiliary operators and the lattice spins to construct a proper fermion representation of the physical 
Hamiltonian. The transform is especially useful for lattice spin Hamiltonians, where two-spin interactions of 
XK-type are either absent or exist only within one-dimensional chains and where the chains are coupled via 
two-spin interactions of Ising type, ring-exchange terms, or more general multi-spin interactions that involve 
an even number of spin operators from each chain. Using the proposed fermionization method we provide a 
simple argument suggesting that a spin Hamiltonian closely-related to the ring-exchange model proposed by 
Paramekanti, Balents, and Fisher [Phys. Rev. B 66, 054526 (2002)] indeed realizes a spin-liquid state. 



Quantum magnetism is one of the richest areas of re- 
search in condensed matter physics. It is also one of the 
most complicated and least understood subjects, primarily be- 
cause most quantum magnets represent strongly correlated 
systems, which do not necessarily relate to a simple non- 
interacting model that could be studied using perturbative 
techniques. The phases and low-energy excitations of quan- 
tum magnets are non-universal and much is determined by 
energetics. Some insight can be obtained using effective field 
theories which represent low-energy physics of various 
low-temperature magnetic phases, but to obtain exact unam- 
biguous results is rarely possible. 

One-dimensional quantum magnets stand out as a notable 
exception to this state of affairs. In many cases, quantum spin 
chains can be solved exactly using various non-perturbative 
methods Perhaps, the simplest such example is a chain of 
half-integer spins with nearest-neighbor interactions, i.e., 

"HxY, Id = -JTjn \^n^n+\ + h.c.j. The nature of the ground 
state of this spin model can be understood with the help of 
Jordan-Wigner transformation, proposed back in 1928 [3]. 
Since our Letter intends to generalize it to some higher- 
dimensional models, let us reiterate the key general ideas, 
which of course are well-known: To solve a quantum spin 
Hamiltonian implies to either calculate the partition function 
or to prove that it is equivalent to a partition function as- 
sociated with a different Hamiltonian, which we understand 
well, preferably a Hamiltonian expressed in terms of cre- 
ation/annihilation operators of some canonical fermions or 
bosons. The main difficulty in accomplishing this task is due 
to the fact that the spin operators, &^ - {crl ± /tr^) /2, are 
neither fermions nor bosons, because they commute on dif- 
ferent sites l^crJJ, (t* j oc §„„ and anti-commute on the same 

site: 1^(7^, (T^j = (here and below Latin indices, a,b - +). 
Clearly, no local transformation can "correct" the anticommu- 
tation relation, but Jordan and Wigner showed that there exists 
a non-local transform, now bearing their names, that accom- 
plishes just that ||3J 



where and /„ are creation/annihilation fermion opera- 
tors, [Z,,/;*,]^ = 6nm- Since = 1, {a^„f = 0, and 

07, - 2(T,^o"^ - 1, the Hamiltonian for the XY-spin chain 
with nearest-neighbor interactions becomes that of a one- 
dimensional Fermi gas in the language of /-fermions The 
Heisenberg spin chain, or a more general XYZ-spin chain 
with nearest-neighbor interactions, take the form of interact- 
ing fermion models, which we recognize as a Luttinger liquid. 
These mappings are possible because the infinite products or 
strings in Eq. ([T) collapse into unity operators in all terms of 
these Hamiltonians. 

It is clear that transform ([T]i itself is independent of the un- 
derlying spin model and it always produces proper fermion 
operators out of the original spins. The problem however 
arises when we use these fermions to describe a generic spin 
Hamiltonian, e.g., a higher-dimensional lattice model [01 or a 
spin chain with longer-range interactions: The string products 
there generally do not disappear and lead to a Hamiltonian, 



which features terms oc /„' exp 
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parameterization- and model-dependent set of sites surround- 
ing a site n. These remaining strings lead in general to vari- 
ous fermion-gauge-theories with constraints. These are com- 
plicated models that we generally do not know how to solve 
and this renders the usefulness of the canonic Jordan-Wigner 
fermionization doubtful. In this Letter, we suggest that by en- 
larging the Hilbert space of the original spin model to include 
"external" auxiliary states, one can fermionize a large class of 
quantum spin models in such a way that both the gauge/string 
factors and the new auxiliary operators disappear from the re- 
sulting fermion Hamiltonian, which therefore takes the famil- 
iar form of an interacting fermion model. Before proceeding 
to specific examples, let us first present the general idea of this 
construction. 

Consider a lattice spin Hamiltonian, "Kj expressed 
in terms of the Pauli matrices, which satisfy the familiar com- 
mutation relations, l^o"", (7^,j - 2i£aj3y&lS„i„- Let us also in- 
clude in the model another Hamiltonian, "H^ [{t/}], which is 
expressed in terms of some other operators, t, that commute 
with all lattice spins and that belong to a representation of an 
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algebra, which may or may not be 5i\{2)k (in all examples be- 
low, we consider spin-one-half systems only). The combined 
Hamiltonian therefore is 

•7/ = +H[{t/}]. (2) 

The main idea here is to use both the original spins and the 
operators of the auxiliary algebra to construct fermion oper- 
ators f„ = F [{(tJI), {t/}] with the goal of achieving the sim- 
plest possible form of the resulting fermion theory. Further- 
more, since the auxiliary algebra may be needed only to "cor- 
rect" the commutation relations, the actual Hamiltonian, "K^, 
in Eq. ^ can be assumed to take the simplest possible form, 
which quite generally implies that "H^ = 0. Hence the parti- 
tion function of the original spin model Zj - Tr exp^-yS"?/,) 
is simply given by Zj - Z/Zj, where Z^ - dj is the dimen- 
sionality of the auxiliary Hilbert space, which is an irrelevant 
constant independent of any parameters of the physical model 
or temperature. On the other hand, a fermionization transform 
(unspecified at this stage) would lead to the Hamiltonian, 

'f(^Hs[{fl}Afn}ATl]]+0. (3) 

We show below that there exists a large class of quantum spin 
models that can be conveniently fermionized with the help 
of auxiliary "external" operators, associated with free stand- 
alone spins [while the Jordan- Wigner transform ([T]i that in- 
volves only "internal" spins is generally not as useful]. 

To provide the simplest explicit example, let us consider 
first a two-leg spin ladder with the following Hamiltonian 

•HxY,2ieg = ^ {[J\\m\n, l)&-{n + 1,0+ h.c] - Ba-"(n, I)] 

«;/=l,2 

+ ^ [7,(n,m)o-^(n,Zi)(T^(m, /a) +h.c.] . (4) 

Eq. (|4]i describes two XF-spin chains in a magnetic field, cou- 
pled via Ising-type interactions. Let us now introduce an aux- 
iliary spin- 1/2 particle, f, with zero Hamiltonian and define 
the following composite operators 

f^(n) = (T^(«, 1) Y\ [o-'(m, 1)] f, and 

f^(n) = <T^(n, 2) Y] 2)] f„ (5) 

where here and below we identify = / and = = 

(/-)'. Since [&~(n,l)f = [r'f = 1, the inverse transform is 

simply a-=^(«, 1 / 2) = fyAn) Um<n {&^{m, 1 / 2)] f.,/^. The 
infinite products in Eq. ^ represent the familiar Jordan- 
Wigner strings ([1]) that give rise to the desired anticommu- 
tation relations within the chains, while the last factor restores 
proper fermion algebra for all operators involved, i.e. 

where ni_2 e Z labels sites, Z1.2 - 1,2 =t, i labels legs, 
and fli_2 = ± distinguishes creation and annihilation opera- 
tors. Furthermore, since d-~{n, I) = 2(t^(«, l)&'{n, I) - 1 - 



2fl(n)fi(n) - l,Vn,l, the fermion representation of (|4|l does 
not involve the auxiliary r-operators (c.f., Refs. llsl^]): 

•WxY.2ieg = - J] {[/||(0/; («)//(« + 1) + h-cj + B/7(n)/;(n)} 
+ [Mn, m) [2f\{n)f^{n) - l] [2fl{m)fy(m) - l] + h.c.} ,(6) 

n.m 

where we have omitted an unimportant constant. We see that 
transform (|5]l maps the spin model onto that of two species 
of fermions with interactions. Let us note here in passing 
that the inclusion of ferromagnetic couplings would translate 
in Eq. (|6]l into an attractive interaction between the fermions, 
which is expected to induce pairing correlations at T ^ 0. In 
particular, the nearest-neighbor intra-chain attraction would 
select an i-wave pairing, while a combination of nearest- 
neighbor and zigzag intra-chain couplings ^ are expected 
to give rise to a p-wave Cooper pairing ITJ for fermions in 
Eq. (|6]l. The zero-temperature ground state of the latter model 
may be related to a one-dimensional topological superconduc- 
tor The spinless version of this phase was in fact one of the 
first models proposed by Kitaev [8] as a realization of a topo- 
logical qubit, which would be associated here with the Jackiw- 
Rebbi/Majorana modes 19] at the ladder boundaries. 

The proposed fermionization scheme can be generalized to 
more complicated Hamiltonians and in particular to higher- 
dimensional lattice spin models. The only requirement is 
that in order for the auxiliary operators to disappear from the 
resulting fermion Hamiltonian, the original spin interactions 
should involve an even number of spin operators from each 
leg. Note also that a "leg" does not necessarily have to be an 
actual row or column of a lattice, these legs may be formed by 
any one-dimensional paths, as long as the collection of these 
paths covers the entire lattice and no two paths intersect. 

Let us now take the next simplest step and consider a 
periodically-translated two-leg spin ladder (|4]i to form a two- 
dimensional square lattice consisting of coupled XY-spin 
chains. We now define the following fermion operators as- 
sociated with the lattice sites, r - in, I): 

fin, D = &Hn, I) W [d^im, I)] (f^) and 

m<n k<l 

fin, D = a-in, I) [crim, I)] (f^) ff , (7) 

where, we introduced a spin- 1/2 operator for each leg [la- 
beled in Eq. O by indices k and /]. The existence of these 
new operators does not disturb the anticommutation relations 
within each chain, because for any product (ni, /)/"'(n2, 
the "new" r-strings square up to an identity operator, while 
the "old" cr-strings enforce the usual constraints. However, 
the existence of the r-operators is crucial to produce the right 
anticommutation relations for operators in difi'erent legs, e.g., 
for /] < h, we obtain: 

[rinuh),r-in2,l2)]^^d-'"inuli) f] mmJi)] (8) 
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FIG. 1: (Color online): Displayed is a square lattice associated with 
the general spin model l|9) and also shown are allowed and forbidden 
cycles on the lattice (depicted as green and red loops correspond- 
ingly). Each loop represents a spin interaction term in the Hamil- 
tonian, which involves a product of L (r*(r)-operators, where L is 
the number of sites in the loop, Cl- Shown here are some exam- 
ples: (i) C2 corresponds to a nearest-neighbor XK-coupling; (ii) C4 is 
associated with ring exchange terms such as in Eq. I ll lb . etc. 



m2<«2 h<k<l2 

Therefore, operators ^ are indeed fermions. Per the same 
arguments as in the two-leg-spin-ladder case (|4]i, we find that 
its two-dimensional realization can be fermionized via (|7]i in 
such a way that no Jordan-Wigner strings, nor r-operators ap- 
pear in the Hamiltonian: We simply can replace &^{n, I) — > 
/"'"(«,/), o--(n,0 ^ f(n,l), &Hn,l) = \lf' {n,l)f{n,l) - l\ 
and change the sign of the XF-interaction term to account for 
&^{n, l)d^{n, I) = -o-+(n, /). 

We now can construct a general spin model in two dimen- 
sions for which the particular transform ^ directly applies. 
Let us reiterate that by "applies," we mean that the result- 
ing fermion Hamiltonian does not contain the r-operators or 
any phase factors arising from Jordan-Wigner strings. Other- 
wise, the transform (|7]) as well as ([1]) always apply in the sense 
that they produce the fermion operators using the spin opera- 
tors available. Let us write down the following rather general 
Hamiltonian associated with a square lattice 



(9) 



ik = h, [{a-;}] + [^"(n, l)J,p{n, l)a^(n + 1,1)] 

n.l 

+ Z Z ■/[Mr)KC]]~[<T''«(r) + h.c. 

CcEven loops |fl(r)=±) reC 

and decipher the meaning of each term in Eq. The 
first term is an arbitrary function of (x'-operators on each 
site. The second term describes a collection of decoupled 
spin chains with arbitrary nearest-neighbor interactions within 
each chain. The last term includes all "allowed" interactions 
that involve two or more different chains (apart from those 
interactions, which may already be present in /z,). The in- 



dices fl(r) - + there label either lowering or raising op- 
erators, which may in principle appear in any combination, 
if no global spin conservation is imposed (no number con- 
servation/gauge invariance for the fermions). Finally, C in 
Eq. (|9]l denotes a closed path, or a one-cycle, on a lattice: 
Cl-{t^\ — > r2 — > r3 — > . . . — > ri], where, r,-, represents a lat- 
tice site and |r,+i - r,| - |r£ - ri| = or 1. The label "Even 
loops" in the sum correspond to the cycles that satisfy the 
following two constraints: (i) Each even cycle has an even 
number of sites present from each leg and (ii) For any site 
r = (n, /) 6 C lying in the path, there exists a r' = («', Z) £ C 
from the same leg, which is either the same site in' - n) or its 
nearest neighbor (n - n' + 1). We note that including the first 
two terms in Eqs. ^ along with the general third term is in 
fact redundant, because the (T~(r)-terms can be associated with 
Cz = {r — > r) cycles (see. Fig. 1) and the nearest-neighbor XF- 
interactions with C2 = {r — » r + e^}. The purpose of all these 
complicated constraints is simple: We want to avoid both r- 
operators in the resulting Hamiltonian [recall that (f'l''-')^" - 1, 
hence the choice of even cycles] and any phase factors that 
could arise from the Jordan-Wigner string (hence, we elimi- 
nate the non-nearest-neighbor interactions within legs |10]). 
A general Hamiltonian (|9]) satisfies all these conditions and 
can be fermionized via the proposed transform (|7]i. Fig. 1 
illustrates the imposed constraints more clearly: The green 
loops there are examples of "allowed" paths, while the red 
loops provide examples of interaction terms, which are not 
allowed. Let us note that there may be a simple physical ex- 
planation why a Hamiltonian with terms of type Cg probably 
can not be fermionized via any variation of our method: The 
loop encircles a site that does not belong to it [(n = 8, / = 6) 
in Fig. 1] and this site may or may not contain a fermion, 
which leads to an ambiguity in the phase factors. Hence, the 
choice of statistics matters in this case, while for Hamilto- 
( 9k it may become largely a question of convenience (c.f.. 

El). 

We can generalize this scheme further by considering 
more complicated long-range intra-chain interactions and var- 
ious lattice structures including three-dimensional and higher- 
dimensional lattices (which too can be conveniently fermion- 
ized using this method as long as the interaction terms respect 
the constraints, i.e., belong to the allowed one- and higher- 
dimensional cycles). However instead of generalizing Hamil- 
tonian (|9]l, let us on the contrary focus on the following spe- 
cific realization of it, which involves very interesting ring- 
exchange terms (associated with loops of type C4 in Fig. 1): 

•HsL^Z {[J\\^^(r)^'(r +e.^ + h.c] - B&-\r)] +ik,+^dm 

reZ- 

where 'H,- J.Yj [o"'(r)(T'(r 4- e^) + h.c.j, c.f. Eq. ©, and 
•7Tf4 = Z o"^(r)o""(r + e.t)o-^(r + e.v + ty)&'(Y + ty) + h. c. 

r 

(11) 

is a ring-exchange term. Our method immediately gives the 
fermionized version of ( fTOl i in the form of an array of Lut- 
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linger liquids coupled via density-density and current-current 
interactions. It is known that such a model hosts a vari- 
ety of unusual phases, including so-called sliding Luttinger 
liquids lll2l4l4ll . This immediately implies that the underly- 
ing spin model too remains a spin liquid in the correspond- 
ing parameter range. Note that Eq. ( fTOl l is closely related 
to the model proposed earlier by Paramekanti, Balents, and 
Fisher y.5] as a realization of a sp in/Bose-liquid. Indeed, 
the rotor Hamiltonian of Ref. Hj, = VZinr-nf - 

r 

K 2 cos (^A^y^r), reduces in the hardcore-boson-limit to a 

spin-(l/2) model of type ( fTOl ) with inter-chain hopping terms 
absent (i.e., J\\ = 0). Since the latter model can be fermion- 
ized as shown, the appearance of a liquid phase predicted in 
Ref. lIlSll is quite natural. Let us also note that model ( fTOl i is 
related to a Bose-metal phase introduced by Motrunich and 
Fisher [16]. The proposed spin/Bose-metal Hamiltonian is of 
the following type "Kbm = Z o-(r)-d-(r')+'H4, where (r, r') 

<r,r'> 

corresponds to nearest-neighbors on a square lattice or a spin- 
ladder 1 3 12]- It is conceivable that the Bose metal may 
represent an isotropic version of the anisotropic spin liquid, 
which appeared in the array of XY-spin chains coupled by ring 
exchanges ( fTOl i. However, the presence of two-dimensional 
XF-couplings in the generic Bose-metal Hamiltonian does not 
allow us to use the fermionization transform O. 

Let us note that Jordan- Wigner-type fermionization with or 
without external operators, when applied to two- and higher- 
dimensional isotropic XY-models would produce a gauge 
theory for the fermions as opposed to a simple interact- 
ing theory, such as in Eq. (|5]|. However, there may exist 
an alternative way to fermionize the X7-model in terms of 
BCS fermions with spin instead of spinless Jordan-Wigner 
fermions. It is known that the partition function of the 
Tgll of interacting fermions with the 



Richardson model 1118 , 
Hamiltonian, "Hr - 2 



Z -/(r,r')c 



r,.!=T.i 

can be mapped onto a combination partition functions, Zr = 
2 (Zspin X Zp_spin), of real spins on singly-occupied sites 
and Anderson pseudo-spins 02011 on paired/empty sites, with 
Zp_spin associated with an XF-model for Anderson pseu- 
dospins, ikxY - T, -/(r, r') a"rO",, + h.c. . Hence, by adding 

r,r' 

to Richardson Hamiltonian non-linear terms that eliminate 
single occupancy, we can suppress the real-spin-sector and 
thereby fermionize the remaining X7-model associated with 
the pseudospins. 

In conclusion, let us mention that the general proposed 
scheme of using "external" operators to build a fermionized 
version of a spin model does not necessarily need to in- 
volve the particular transform (|7]i, which should be viewed as 
merely an example. Furthermore, there may exist multiple dif- 
ferent parameterizations or numberings of both the "external" 
operators and "internal" Jordan- Wigner-like strings that give 
rise to the same Hamiltonian and associated partition function 
(modulo an overall constant). To provide a final example in 



this context, let us consider a three-leg spin ladder, forming a 
triangular lattice. To fermionize this model, we in addition to 
transform (|7]i, can use, e.g., the following fermion operators, 
fj^(n) = o-^(n, l)nm<« [d-~(m, l)]f', where I ^ 1,2,3 = x,y,z 
is the leg index. Hence, we can fermionize this three-leg lad- 
der using just one auxiliary spin, f . However, whether or not 
this particular parameterization is useful depends on a Hamil- 
tonian of interest. If only nearest-leg ring exchanges are in- 
cluded, this transform leads to a proper fermion Hamiltonian, 
which does not involve the external spin. However, the in- 
clusion of center-of-mass conserving ring-exchanges that in- 
volve all three legs (such as present in the models discussed 
in Refs. ll2lll22ll ) would lead to the appearance of f*-"gauge"- 
factors in the fermion Hamiltonian. 
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